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An analysis of CoRoT multicolour photometry of exoplanets 
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ABSTRACT 

We analysed the chromatic data of the planetary transits observed with CoRoTXo supply ho- 
mogeneous time series in each of the CoRoT colours. In a first step, we cleaned the chromatic 
light curves from the contamination of nearby stars and removed outliers and trends caused by 
anything different from the planetary transits. Then, we simultaneously fitted the chromatic 
transits of each planet, obtaining a common solution for the orbital parameters i, to and a/R s , 
with a particular care in the fitting for different limb-darkening parameters. The planet-to-star 
radius ratios in the CoRoT colours are compatible when considering one planet at a time, but 
the ensemble of low-mass planets seems to show a peculiar behaviour of the radius ratios. 

Key words: methods: data analysis - techniques: photometric - planetary systems. 



1 INTRODUCTION 

The discovery of new planets is recurrent news in the exploitation 
of ground-based surveys (e.g., Mayor et al. 2011) and space mis- 
sions (e.g., Borucki et al. 201 1). We count more than 700 confirmed 
detections and a lot of candidates are in the waiting list. Therefore, 
it is quite natural to direct our attention to techniques that can help 
in disentangling the uncertain cases and/or supply more informa- 
tion on the planet's structure. If we are able to point out any kind of 
new feature, this can help in planning future investigations and/or 
theoretical analysis. 

Kepler and CoRoT (convection, rotation and planetary transits; 
Baglin et al. 2006) space missions make use of the photometric 
method of transits to detect exoplanets. Theory predicts that the 
depth of the transit of a planetary body in front of its parent star 
should not depend on the wavelength, except for the effect of the 
stellar limb darkening. However this is true only at first order: the 
apparent radius of the planet, and therefore its transit depth, does 
depend on the wavelength; the light from the parent star filters 
through the planet's atmosphere, so that the atmospheric signatures 
remain imprinted in the transit depths at different wavelengths, cre- 
ating a transmission spectrum. The weight of this effect depends 
on the atmospheric scale height. After the first Hubble Space Tele- 
scope positive result in 2002, the detection of spectral features in 
the exoplanet atmospheres over a broad wavelength range is now a 
firmly established field (Seager & Deming 2010). 

In the case of the CoRoT mission, a chromatic device was 
mounted in the focal block with the aim of discriminating 'nearly 
achromatic planetary transits from coloured stellar fluctuations' 



(Barge et al. 2006). The light passes through a prism before en- 
tering in the exoplanet channel: to obtain chromatic light curves, 
the mask is divided into three submasks along the dispersion direc- 
tion. The mask (and so the relative submasks) for each target has 
a shape optimized for the magnitude and location of nearby fainter 
stars. Such a shape is chosen between a set of templates in order 
to maximize the signal-to-noise ratio (S/N). The chromatic device 
plays an important role in the selection of the most promising tar- 
gets, allowing the rejection of several false positives (Carpano et al. 
2009), and the confirmation of the planetary nature of objects like 
CoRoT-7b (Leger et al. 2009) and CoRoT-8b (Borde et al. 2010). 

The specificity of each mask constitutes a great difficulty in ex- 
tracting physical parameters from the coloured light curves, but the 
efforts involved could be rewarding. Indeed, the detection of small 
differences in the amplitudes and in the phases of the chromatic 
light curves of the high-amplitude 6 Set star CoRoT 101 155310 al- 
lowed us to discriminate between the radial and non-radial nature of 
the pulsation mode (Poretti et al. 2011). Therefore, we considered 
it noteworthy to undertake a study of the chromatic light curves of 
the first planets discovered by CoRoT to point out the useful hints 
that this unique photometric device can supply for further investi- 
gations. 



2 DATA ANALYSIS 

At the time of our analysis, CoRoT has discovered 23 planets, of 
which only the first 15 have public data. Among these 15 planets, 
only 10 belong to the subset for which chromatic data are available: 
our analysis focused on these targets (Table[T). 



* The CoRoT space mission was developed and is operated by the French 
space agency CNES, with participation of ESA's RSSD and Science Pro- 
grammes, Austria, Belgium, Brazil, Germany, and Spain, 
t E-mail: francesco.borsa@brera.inaf.it 



2.1 Removing contamination 

The point spread function (PSF) of the CoRoT optics is highly de- 
focused, to improve the photometric precision and to reduce the 
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Table 1. Estimates of the contamination factors in each mask (white) and submask (red, green, blue) for the 10 CoRoT planets analysed, and characteristics of 
the parent stars. Nearby stars have been considered up to a 20 arcsec distance from the source. 



Planet CoRoT colour contamination Parent V Contaminants Discovery 

White Red Green Blue Star mag Num. V mag range paper 

(per cent) (per cent) (per cent) (per cent) 



CoRoT-lb 
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31 
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36 +0 
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1 1 


16+° 07 

"■ lo -0.06 


CoRoT-2b 
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00 
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02 
02 
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CoRoT-8b 
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CoRoT-9b 
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CoRoT-llb 
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65 
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i 7+0.99 
-0.82 
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08+° 


37 
23 


1 Q0+ 044 

l y -0.34 


CoRoT-3b 
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02+° 


97 
84 
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1 9+0.45 
1? -0.41 


11 


.14+ 


.96 
.84 


17 60+ 3 49 

"'""-2.74 



sensitivity to satellite jitter. Such defocusing increases the proba- 
bility that the PSF of the star harbouring the planet overlaps the 
PSFs of nearby stars. None of the N2 coloured light curves is cor- 
rected for the fluxes of stars entering the target star mask. There is 
only an estimate of the contamination factors made with a generic 
mask before the run, but only for the white light curve. An analysis 
of the contamination effects in the CoRoT colours was performed in 
the case of CoRoT-8b only (Borde et al. 2010). Since the quantita- 
tive evaluation of the contamination effect is necessary to compare 
transit depths in different colours, we applied a similar procedure 
by using the data available in the EXODAT database (Deleuil et al. 
2009). 

A preliminary step was the evaluation of the projected image on 
the CCD of the contaminant stars only, without the main target. We 
calculated the inflight PSF for each target and then we applied it 
to the nearest six contaminant stars (i.e., those listed in EXODAT), 
scaling the calculated PSF to match their V magnitudes and B - V 
colour indices. If the six EXODAT contaminants were all closer 
than 20 arcsec to the target, we completed the list with stars taken 
from the USNO-A2 catalogue, in order to be sure to include all 
the potential contaminating sources. Fig.Q]shows the histogram of 
the magnitude differences AV between the 67 contaminants and the 
stars harbouring planets. In two cases only we have AV < 2.0 mag, 
related to contaminants far away from the main targets, i.e. 17.2 
and 18.7 arcsec. This is not surprising, since the masks are actually 
chosen to minimize the contamination effect. However, although 
very small, this effect cannot be neglected for our purposes. 

The above procedure supplied us a clear starting point, i.e. it al- 
lowed us to evaluate how many pixels of the target mask were con- 
taminated by the flux from other stars, on the basis of their position, 
magnitude and colour. To proceed in the photometric analysis, we 
had to take into account not only the different colours of the stars 
but also the impossibility of having accurate photometric transfor- 
mations between the V and B - V values and the spectral response 
of the specific mask+star combination. 

To tackle this severe, intrinsic problem of the CoRoT photometry, 
we divided the calculated PSF of each contaminant into two parts, 
one bluer and one redder. We scaled the bluer part as the B mag and 
the redder part as the V mag with respect to the main target. 

Then, we integrated the flux of the contaminants in the main 
mask and in the three submasks, comparing it with the flux pro- 
duced by the main target. To estimate the contamination factor and 
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Figure 1. Histogram representing the V magnitude differences between the 
main targets and the contaminant stars. Bins refer to intervals of 0.5 mag. 



its error bar we repeated the computations by running the proce- 
dure 20 000 times for each star. We left the parameters vary in wide 
ranges: the B and V magnitudes by ±0.3 mag (and hence the B—V 
colour by +0.6 mag) and the center of the image of each star in a 
3x3 pixel 2 box. The choice of these wide ranges was done for sev- 
eral reasons. EXODAT does not provide information about uncer- 
tainties on both magnitude and position of the contaminants, which 
are very faint stars. Moreover, the large error bar on the colour in- 
formation is a way to take into account the impossibility of using 
the exact CoRoT spectral response on contaminants. The magni- 
tude range also accounts for possible light variability. At the end of 
the process we chose the median values as the best estimate of the 
contamination factors and their error as the 68 per cent confidence 
interval (Table [T}. The contamination factors for CoRoT-8b are in 
excellent agreement with those estimated by Borde et al. (2010), 
i.e. 0.9, 0.7, 0.2, and 2.4 per cent for the white, red, green and 
blue CoRoT colours (hereafter W c , Rc, Gc, B c ), respectively. This 
agreement strengthened our confidence in the adopted procedure. 

The green submask is always composed of only one column of 
pixels: it acts as a separator between the blue and red spectra. As 
a consequence, the green flux is low compared to the others, and 
suffers from a higher dispersion. Hence we decided to disregard 
the Gc light curve, concentrating only on the Rc and Be channels. 
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Table 2. List of the CoRoT planets analysed, with relative periods, 
ephemerides and masses taken from the discovery papers listed in Table[T] 



Planet 


Period (d) 


t - 2450000 


Mass (Mj) 


CoRoT-lb 


1.5089557 ±0.0000064 


4159.4532 


1.03 


CoRoT-2b 


1.74299641 ± 0.0000017 


4706.4041 


3.31 


CoRoT-4b 


9.20205 ± 0.00037 


4141.36416 


0.72 


CoRoT-5b 


4.0378962 ± 0.0000019 


4400.19885 


0.467 


CoRoT-6b 


8.886593 ± 0.000004 


4702.2556 


2.96 


CoRoT-7b 


0.853585 ± 0.000024 


4398.0767 


0.015 


CoRoT-8b 


6.21229 ± 0.00003 


4238.9743 


0.22 


CoRoT-9b 


95.2738 ±0.0014 


4603.3447 


0.84 


CoRoT-llb 


2.99433 ±0.000011 


4597.679 


2.33 


CoRoT-3b 


4.25695 ± 0.000005 


4283.1383 


21.66 



2.2 Cleaning and detrending light curves 

After having corrected each light curve with the contamination fac- 
tor estimates, we analysed all the light curves separately, because 
trends and outliers are different in each colour. 

We did not consider observations acquired when the satellite 
crossed the South Atlantic Anomaly and those flagged as inaccu- 
rate by the pipeline. After that, we analysed the CoRoT time series 
to study each planetary transit. The light curves are not linear: they 
show trends due to the ageing of detectors (Auvergne et al. 2009) 
and jumps due to hot pixels hit by cosmic rays. After the removal 
of obvious outliers, we calculated the average flux difference be- 
tween two consecutive points. We deleted the points which were 
more than three times this quantity both from the previous and the 
next points. 

Several algorithms have been proposed to detrend CoRoT light 
curves (e.g. Mazeh et al. 2009). We applied a specific approach 
for our purpose, i.e. having well-cleaned curves around the ob- 
served transits. We selected two intervals in the out-of-transit part 
of the light curves, before and after the central time of each transit 
T u where the light curve is expected to be flat, but actually is af- 
fected by stellar activity and instrumental trends (Fig. [2] top panel). 
As a rule, we considered the intervals [Tj - 2.0 D, T f - 0.7 D] 
and [Ti + 0JD, T t + 2.0 D], where D is the full duration of the 
transit. The data in these two intervals were fitted by a second- 
order polynomial. Then, we considered all the data in the interval 
[Tj - 2.0 D, Tj + 2.0 D], and we divided them by the second-order 
fit, thus removing the trend and normalizing the flux units to the 
out-of-transit level (Fig. [2] bottom panel). In such a way, we also 
corrected for stellar activity. The calculated errors on all the de- 
trending parameters have been propagated through the entire anal- 
ysis; the errors on the orbital periods cited in the literature have 
also been taken into account when phase-folding the transits. We 
then superposed all the transits of each light curve using the pe- 
riods determined in the discovery papers (Table |2]l to obtain the 
phase-folded transits in each colour. 

2.3 Limb darkening and transit fitting 

The treatment of limb darkening in the fitting of planetary tran- 
sits has long been analysed. Strong efforts were made to calcu- 
late the theoretical limb-darkening parameters, depending on the 
characteristics of the parent star and on the bandpass of the mea- 
surements (e.g. Claret 1998). With the coming of the CoRoT and 
Kepler space-based observatories, new effort was necessary to cal- 
culate new parameters for their ranges of sensitivity (Sing 2010; 




2S05 25 1 

time [HJD-2451545] 




Figure 2. Example of how our algorithm works. Top: raw white light curve 
of CoRoT-2b, in grey are evidenced the data points we considered for the 
second order fit before and after each transit. Bottom: the final phase-folded 
transit of CoRoT-2b after the correction algorithm. 



Claret & Bloemen 2011). However, no estimation has ever been 
made for the CoRoT colours. The CoRoT colours in fact are not 
bandpasses, because they are split by a low-dispersion device, and 
so the exact range of wavelengths slightly varies for each star, de- 
pending on its colour and position on the CCD; this lack of a precise 
photometric system is preventing the theoretical estimation of the 
limb-darkening coefficients. Howarth (2011) stresses how difficult 
it is to use standard and/or theoretical limb-darkening parameters 
when studying exoplanetary systems. We then decided to assume 
the general quadratic law for the limb-darkening effect only and to 
compute the /Ji and /j 2 coefficients as free parameters in each case. 

The phase-folded coloured light curves were fitted with the tap 
package (Gazak et al. 201 1), which uses the Mandel & Agol (2002) 
model for the transit shape. We performed a Markov Chain Monte 
Carlo analysis of 1000 000 steps for each planet, fitting simulta- 
neously the 3 light curves (Wc, Rc, Be) keeping the orbital period 
fixed. The inclination angle i, the distance between the planet and 
the star a (in units of the stellar radius R s ) and the transit timing 
to were calculated considering all the colours simultaneously. On 
the other hand, the planet-to-star radius ratios R p /R s (the square 
root of the transit depth) and the limb-darkening coefficients \i\ and 
H2 were calculated for each colour (Table [3} . The limb-darkening 
parameters tabulated by Sing (2010) for the CoRoT range of sensi- 
tivity are within the error bars of our values for the Wc light curves 
(Tabled, thus corroborating our methodological approach. The re- 
ported error bars are calculated by propagating the errors on the 
parameters (the long-term trend in different colours, the orbital pe- 
riods and the contamination effects) we used to extract the shape of 
the transit from the light curves. 

Regarding the orbital parameters, all orbits were considered to be 
circular, except for CoRoT-5b (e = 0.07 + 0.06) and CoRoT-9b (e = 
0.07 ± 0.05), for which eccentricities were left as free parameters 
but constrained around the values found in literature. The resulting 
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Table 3. Fit results for the phase-folded transits of the planets analysed, with planet-to-star radius ratios and limb-darkening coefficients (linear and quadratic) 
for the Rc and Be colours. The orbital periods are taken from the discovery papers. 



Planet 


a/R s 


i 


RpIRs 






J"2 








o 


R C 


B C 


R C 


B C 


R C 


B C 


CoRoT-lb 


4.72 ±0.17 


83.8 ± 1.0 


0.1440 ±0.0028 


0.1471 ±0.0030 


0.58 ±0.31 


0.34 ± 0.27 


-0.16 ±0.45 


0.38 ± 0.38 


CoRoT-2b 


6.76 ± 0.09 


88.06 ± 0.72 


0.1655 ±0.0015 


0.1665 ±0.0021 


0.38 ±0.10 


0.55 ± 0.20 


0.03 ± 0.21 


0.06 ± 0.33 


CoRoT-4b 


16.60 ± 1.05 


89.01 ±0.75 


0.1067 ±0.0025 


0.1052 ± 0.0041 


0.40 ± 0.29 


0.55 ±0.31 


-0.01 ±0.41 


-0.03 ± 0.44 


CoRoT-5b 


9.44 ± 0.77 


85.55 ± 0.72 


0.1 144 ±0.0030 


0.1 104 ±0.0050 


0.45 ± 0.30 


0.54 ± 0.35 


0.00 ± 0.39 


0.08 ± 0.47 


CoRoT-6b 


17.96 ±0.89 


89.04 ±0.61 


0.1152 ±0.0024 


0.1174 ±0.0031 


0.29 ± 0.23 


0.51 ±0.29 


0.20 ± 0.41 


0.00 ± 0.45 


CoRoT-7b 


5.40 ± 1.10 


84.8 ±4.5 


0.0182 ±0.0015 


0.0158 ± 0.0088 


0.63 ± 0.29 


0.45 ± 0.35 


0.01 ±0.39 


0.00 ± 0.38 


CoRoT-8b 


18.3 ±2.7 


88.56 ± 1.03 


0.0811 ±0.0034 


0.0755 ± 0.0072 


0.58 ± 0.30 


0.50 ± 0.32 


-0.07 ±0.41 


0.04 ± 0.39 


CoRoT-9b 


96.5 ±5.0 


89.84 ±0.12 


0.1 144 ±0.0026 


0.1202 ± 0.0042 


0.42 ± 0.32 


0.42 ± 0.34 


-0.11 ±0.48 


0.07 ± 0.44 


CoRoT-llb 


6.71 ±0.37 


82.86 ± 0.62 


0.1009 ±0.0034 


0.1031 ±0.0024 


0.44 ± 0.34 


0.52 ± 0.33 


0.05 ± 0.41 


-0.04 ± 0.38 


CoRoT-3b 


8.3 ± 1.3 


86.5 ± 2.6 


0.0686 ± 0.0030 


0.0603 ± 0.0080 


0.52 ±0.35 


0.52 ± 0.34 


-0.15 ±0.49 


0.00 ± 0.39 



Table 4. Comparison between fitted limb-darkening coefficients on the We 
light curves and those estimated by Sing (2010) for the CoRoT range of 
sensitivity. 



Planet 



/'i 



V2 



W c 



Sing 



W c 



Sing 



CoRoT-lb 


0.51 


+ 


0.27 


0.37 


-0.03 


± 


0.40 


0.27 


CoRoT-2b 


0.43 


+ 


0.10 


0.46 


0.02 


± 


0.20 


0.22 


CoRoT-4b 


0.43 


± 


0.36 


0.37 


0.04 


± 


0.52 


0.28 


CoRoT-5b 


0.37 


+ 


0.31 


0.37 


0.11 


± 


0.38 


0.27 


CoRoT-6b 


0.36 


+ 


0.26 


0.38 


0.09 


± 


0.43 


0.27 


CoRoT-7b 


0.48 


+ 


0.35 


0.53 


-0.02 


± 


0.49 


0.17 


CoRoT-8b 


0.52 


+ 


0.31 


0.59 


-0.03 


± 


0.40 


0.12 


CoRoT-9b 


0.54 


+ 


0.25 


0.46 


0.17 


± 


0.37 


0.22 


CoRoT-llb 


0.46 


+ 


0.34 


0.34 


0.02 


+ 


0.40 


0.29 


CoRoT-3b 


0.44 


+ 


0.34 


0.32 


0.00 


+ 


0.43 


0.31 



values of the inclination angle and of a/R s are all within the error 
bars of those given in the discovery papers. 

We also investigated the possibility that the decision to fit si- 
multaneously the limb-darkening parameters could have affected 
the R p /R s determinations. We repeated the same fitting procedure 
keeping fixed the linear /Ji parameters to the median values (0.44 
for Re, 0.51 for Be)- We obtained the same values for the radius 
ratios, with negligible differences with respect to the values listed 
in Tabled 



2.4 The case of CoRoT-3b 

With its mass of 21.66 Mj (Deleuil et al. 2008), CoRoT-3b was 
defined as an inhabitant of the 'brown dwarf desert' . It could be a 
low-mass brown dwarf or the member of a new class of 'superplan- 
ets', because it is on the edge between planets and stars. Therefore, 
we did not consider it as a bona-fide member of our class of exo- 
planets, but the presence of two close stars 2.9 and 4.9 mag fainter 
at distances of 5.6 and 5.3 arcsec, respectively (Deleuil et al. 2008), 
makes its analysis very challenging for the relevance of the contam- 
ination (Table [TJ. 

The light curves of the original chromatic data, without apply- 
ing any correction for the contaminant stars, show a clear differ- 
ence of 0.002 in relative flux in the depths of the transits (Fig. [3] 
upper panel). After applying our method to remove the contamina- 
tion using the photometry listed in EXODAT, this amount is halved 



.002 



.000 



0.998 



0.996 



0.994 



i. - 



J- 



2695.3 2695.4 2695.5 2695.6 2695.7 2695.8 2695.9 
time [HJD-2451 545] 



0.998 



0.996 



0.994 



■ ■ <■ * 



2695.3 2695.4 2695.5 2695.6 2695.7 2695.8 2695.9 
time [HJD-2451545] 

Figure 3. Transits of CoRoT-3b in the Re (triangles) and Be (squares) 
colours, before (top) and after (bottom) removing the contamination of 
nearby stars. Points refer to a binning of 0.003 d. 



(Fig. [3] lower panel), confirming that our procedure works in the 
right direction. The resulting difference between the radii in Re 
and Be is still large (~12 per cent, see Table [3j, but the scattered 
B c light curve and its quite large error bar makes the values of the 
radii compatible in this extremely contaminated case. In Fig. [3] a 
small 'bump' is also notable at the bottom of the B c transit. Due 
to the fact that the orbital period of CoRoT-3b is compatible with a 
synchronization with the stellar rotation (Mazeh & Faigler 2010), 
we cannot exclude it to be an inhomogeneity on the surface of the 
F3V parent star. We also should not forget the specificity of the 
nature of the CoRoT-3 system. For the sake of completeness, we 
calculated that a missed bright contaminating star with V=16.0 and 
B- V=+0.3 could equalize the depths of the transits when affecting 
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the B c contamination up to 35 per cent. Undetected or undetectable 
contaminant stars, especially in the crowded field of the Galactic 
plane observed by CoRoT, constitute a possible error source we 
cannot quantify. However, at the moment, there is no evidence of a 
missed contaminating star around CoRoT-3 (Deleuil et al. 2008). 



3 DISCUSSION 

The planet-to-star radius ratios in the Rc and Be colours (Table[3} 
provide the expected evidence that the transits are mostly achro- 
matic. In several cases, the ratio values are far from equality, but the 
error bars recover it in all cases except for CoRoT-9b. Therefore, 
we looked at the whole set of planet-to-star radius ratios to verify if 
the ensemble analysis could provide some hints of departures from 
equality. To do that we used the distance from the parent star, the 
mass of the planet and the effective temperature of the planet. We 
obtained the plot shown in Fig. [4] when sorting the planet-to-star 
radius ratios in function of the planet's mass. We can see that all 
the planets with a B c radius smaller than the R c radius are confined 
in the M < 0.8 Mj region. We checked if such behaviour could 
depend from the contamination level, but the plot did not change 
by varying the contamination factors up to 30 per cent of the value 
calculated by using the procedure described in Section 2.1. This 
could mean that the Rc radius of the low-mass planets is larger than 
the B c radius, while for massive planets the radii are almost equal 
(Fig- ID- The statistics is very limited at the moment and the fact 
that all the planets with M < 0.8 Mj are below the equality line 
could be fortuitous, also considering that they are those with the 
bigger error bars. For the sake of completeness, we verified if the 
behaviour sketched by Figs.[4]and[5]can reflect some physical rea- 
son. Tessenyi et al. (2012) stated that the height of a planet's atmo- 
sphere can be estimated to be Az = nH, where typically n ~ 5 and 
H is the atmospheric scaleheight. Since H= ^R, the smaller grav- 
ity of low-mass planets results in a larger atmospheric scaleheight, 
and this could enhance the atmospheric features in the photomet- 
ric measurements at different wavelengths. We also verified that 
the observed transit depth differences of some hundreds of ppm are 
consistent with the effects of planetary atmospheres. In the specific 
case of the CoRoT planets, we obtained values of Az in the range 
0.5-5 per cent of the radius of the planet (even more for CoRoT- 
7b). These values of the atmospheric height correspond to upper 
limits on transit depth differences due to an atmospheric absorption 
of 100-1300 ppm. 

We searched for other transiting planets with transit depth mea- 
surements in multiple optical bandpasses. The characterization of 
the planetary atmosphere of HD 189733b (M = 1.138 Mj) sup- 
plied evidences of the presence of high-altitude atmospheric haze 
detected in the visible (Pont et at. 2008; Lecavelier Des Etangs et 
al. 2008; Sing et al. 2011) and probably extended to the infrared 
(Gibson et al. 2012). The transmission spectrum is in good agree- 
ment with that of Rayleigh scattering and the planetary radius re- 
sults to be slightly smaller at long wavelengths, in agreement with 
the results obtained from CoRoT photometry (Fig. |5j- The mass 
of HD 209458b (M = 0.69 Mj) is too close to the value 0.8 My 
to provide a useful check on the dependence of the planet-to-star 
radius ratios on the mass, but we note that Knutson et al. (2007) 
obtained different radius values in function of wavelength in the 



k is the Boltzmann constant, T the equilibrium temperature of the planet, 
/i the mean molecular mass of the atmosphere and g the gravity acceleration. 



0.1 5 


8 

1 b 






a 

2b 


0.10 


□ 

^ A9b 

A 1 

□ 

o D 




a 

1 1 b 


8 

6b 


0.05 










0.00 


7b 








1 


? 




3 



Mass [Mj] 

Figure 4. Planet-to-star radius ratios of the CoRoT planets analysed in the 
Rc (triangles) and Be (squares) colours. Error bars are of the order of the 
point dimension (except for CoRoT-7b which are about twice the point di- 
mension). 
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Figure 5. Radius of the planet in the Be colour with respect to the radius 
in the Rc colour. Blue circles refer to our analysis of CoRoT planets. The 
huge eiTor bars of CoRoT-7b, not even completely inside the plot, are due 
to the very small depths of the transit with respect to the S/N, especially 
in the Be colour. Red triangles refer to literature data (cited in Section 3) 
taken in different optical bandpasses: Rbiue refers to the measurement with 
the shorter wavelength, and R ret j to the measurement with the longer one. 



full range from 293 to 1019 nm. Theoretical transmission spectra 
support such transit depth variations (Gillon et al, 2012, see their 
fig. 7). With a mass of only 0.081 Mj, HAT-P-1 lb could provide a 
significant test. Deming et al. (201 1) pointed out that the radius ra- 
tio is 6 per cent larger in the J band than in the Kepler photometry. 
They are not able to provide an explanation for such a large discrep- 
ancy and suggest a variable contribution of the water absorption in 
the Earth's atmosphere during the observations. 
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4 CONCLUSIONS 

We provided a methodological approach to extract the chromatic 
information contained in the light curves of 10 transiting planets 
discovered by CoRoT. We satisfactorily removed the contaminat- 
ing light coming from poorly known stars inside the pre-defined 
aperture masks. This step is of extreme importance for a coherent 
analysis of transits, as third light can affect in different percentages 
the chromatic light curves and thus bring to incongruous compar- 
isons between the different colours. When comparing transits of 
the same planet observed with different instruments (or at different 
wavelengths), an analysis of the contamination effect is fundamen- 
tal, since third light could lead to different and erroneous transit 
depth (and planetary radius) estimations. 

Then we detrended and cleaned the data producing a homoge- 
nous set of light curves in each of the CoRoT colours. The anal- 
ysis of the CoRoT chromatic light curves confirms that all the 
transit depths for the planets studied in two different colours are 
nearly equal. However, we could put in evidence an interesting be- 
haviour of the planet-to-star radius ratios of the transits in Be and 
R c colours when they are plotted versus the mass of the planets. We 
suggest the possibility of observing departures from transit depth 
equality for planets having M < 0.8 Mj. This hypothesis is still 
rather speculative at the moment, but is not groundless, roughly re- 
flecting the relevance of the atmospheric scale height in low-mass 
planets. The inclusion of the Na D lines in the range of the Rc 
colour (Leger et al. 2009) could also play a role when considering 
that Sing et al. (2008a, 2008b) report on the very strong absorption 
of the Na D lines in the high atmosphere of HD 209458b, around 
589 nm. Of course, the results would have been more compelling if 
the CoRoT photometric system was a standard or a well-calibrated 
one. 

The careful evaluation of the chromatic effect could result in 
a fertile exercise in the case of low-mass planets. High-precision 
multicolour photometry of transits could provide an efficient tool 
to investigate the upper atmosphere of exoplanets, taking into ac- 
count the intrinsic difficulties inherent in accessing these thin layers 
in other worlds' atmospheres. 
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